Deuterons and space-momentum correlations in high energy nuclear collisions 
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Using a microscopic transport model together with a co- 
alescence after-burner, we study the formation of deuterons 
in Au + Au central collisions at y^i = 200 AGeV. It is 
found that the deuteron transverse momentum distributions 
are strongly affected by the nucleon space-momentum corre- 
lations, at the moment of freeze-out, which are mostly deter- 
mined by the number of rescatterings. This feature is useful 
for studying collision dynamics at ultrarelativistic energies. 



Relativistic heavy ion collisions offer the unique oppor- 
tunity to study hot and dense matter under controlled 
laboratory conditions (for recent reviews, see 
However, the particle momentum distributions do not di- 
rectly reveal the properties of the initial dense state of the 
system created in these collisions. The system undergoes 
longitudinal and transverse expansion which greatly af- 
fects the particle momentum distributions. Therefore the 
observation that for a given colliding system and within 
the same kinematic region, the slope parameter^ T de- 
pends on the particle mass Q is of great interest and can 
be exploited to extract information about the reaction 
dynamics prior to the freeze-out stage. It has been found 
that the higher the particle mass, the larger the slope 
parameter is. This mass dependence is the strongest for 
the heaviest systems (Pb + Pb), and vanishes altogether 
in p+p collisions at similar energies l|. 

A related observation is that the size parameters Rt 
and Rhong {Rhong and Rt are the size parameters in 
the beam direction and perpendicular to the beam direc- 
tion, respectively) , which are extracted from two-particle 
correlation measurements, depend on the transverse mo- 
mentum of the pair. The higher the momentum, the 



1 The slope parameter is extracted from the fitting of the 
transverse mass distribution — — f N with an exponential 

mT ami 

function A ■ exp(—rriT/T). A is a normalization parameter. 



smaller the size parameter |g,|7|. This dependence, like 
the slope parameter dependence, is the strongest (in 
RLong) in the heaviest systems, although it is also ob- 
served in elementary collisions S. 

In any heavy-ion collision, the space-time freeze-out 
distribution and its dependence on the particle momen- 
tum are determined by the underlying dynamics. The 
above observations are usually interpreted by considering 
nuclear fluid dynamics (NFD) type collective flow, which 
clearly leads to space-momentum correlations P-pT|. 
The NFD is not the only model for the interpretation 
of such correlation. Another example might be string 
fragmentation p3fr . 

Further insight into transverse collective flow can be 
gained from studying the mean transverse momentum, 
(jot), of different particles like pions, kaons, and protons: 
specifically, analyzing the mean pt variation with particle 
mass. Unfortunately, the full space-momentum structure 
of the collision cannot be extracted from single-particle 
momentum spectra alone To shed more light on 

this issue, we propose to utilize deuteron distributions to 
investigate nucleon freeze-out properties. In this letter 
we use the deuteron transverse momentum distribution 
as well as the ratio of the proton distribution to that 
of the deuteron |l5) , to extract information on space- 
momentum correlations and flow at RHIC energies. In 
this study, we focus on central (impact parameter b < 
3.0 fm) Au + Au collisions at ^/s = 200 AGeV. 

For our investigation we employ the Relativistic Quan- 
tum Molecular Dynamics model (RQMD) The 
model is well established and has been used success- 
fully to describe many observables measured at AGS and 
SPS bombarding energies over a wide range of projectile- 
target combinations. RQMD jl6) is a semi-classical mi- 
croscopic approach which combines classical propagation 
with stochastic interactions. Color strings and hadronic 
resonances can be excited in elementary collisions. Their 
fragmentation and decay lead to production of particles. 
Overlapping strings do not fragment independently from 
each other but form 'ropes', chromoelectric flux-tubes 
whose sources are charge states in higher dimensional 
representations of color SU(3). RQMD is a full transport 
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theoretical approach to reactions between nuclei (and el- 
ementary hadrons) starting from the initial state before 
overlap to the final state after the strong interactions 
have ceased (freeze-out). The model does not include 
light cluster productions. Therefore an after-burner, de- 
scribed in |lq-|2lf|, is used for deuteron yields calculations. 
More details of the coalescence type calculations can be 
found in (l^j2l|,|2^] and references therein. 

The deuteron binding energy (^2.2 MeV) [^3j is small 
compared to the characteristic freeze-out temperature 
(~ 140 MeV) of ultra-relativistic heavy-ion collisions we 
are studying. Hence, deuterons cannot survive rescatter- 
ing. Since many rescatterings occur within the hot and 
dense reaction phase, the only deuterons that survive and 
escape are those formed near the freeze-out stage, either 
on the surface of the fireball or at a later time when the 
environment is dilute. 

The particle phase-space distribution at freeze-out re- 
flects physics at earlier stage of the collision. Similar to 
the two-particle correlation functions, the probabilities of 
bound state (deuterons and heavier clusters) formation 
are determined by this distribution pjj . 

The sensitivity of the two particle correlation measure- 
ment to the source size decreases when the source size 
(and/or duration of the emission) becomes large. The 
analysis of the deuteron yield relative to the proton yield 
has no such loss in sensitivity although the value of the 
deuteron yield decreases due to larger distance between 
the two nucleons. 

Assuming a Gaussian form for the nucleon source 
p6|p7[ , a size parameter R g can be extracted from the 
single particle distributions of protons and deuterons [g8| : 

,jp d a N p s 2 

r 9 = -(V^nc) , (o.i) 

where m p ,md are the proton and the deuteron masses, 
respectively. The invariant distribution is Eid 3 Ni/d 3 p 
with (i = proton, deuteron). The above equation as- 
sumes that the deuteron energy is the sum of the pro- 
ton and the neutron energies and that there is no space- 
momentum correlation in particle distributions at freeze- 
out. Then the space-momentum correlation can be stud- 
ied by inspecting the R g as a function of the transverse 
momentum px- 

The Gaussian size parameter R g as a function of 
the nucleon transverse mass mx is shown in figure 1. 
Here, the filled circles represent the results from origi- 
nal(default) RQMD events. The open symbols represent 
events with altered space-momentum correlation. The 
squares represent the so called aligned case, where for 
each nucleon the space vector fx is aligned with the 
transverse momentum vector fix. The triangles repre- 
sent the case where the angle between fix and fx has 
been randomized. Note that in the aligned and random 
cases only the relative orientation of rr to fix is mod- 
ified: momentum distributions and projections onto ei- 
ther rx or pt are not touched. In the randomized case, 



the amplitudes of vectors \fj<\ and \fix\ are still corre- 
lated. To remove such correlation, the vectors fx and fix 
were scrambled (open circles in Fig.l). After the opera- 
tion, the correlations among fx and fix are removed and, 
as expected, the distribution is almost flat as a function 
of mx- To guide the eye, the solid line represents the 
function 9.75 • (my) 1 / 2 . 

While no large differences are observed between the 
normal and the aligned results, a dramatic effect is evi- 
dent between the randomized and normal cases. This im- 
plies that nucleon momenta are already largely aligned in 
the real events. The calculated mean cosine of the angle 
between transverse space and momentum vector is about 
(cos(0)) ~ 0.9 at mid-rapidity. As one can see in the fig- 
ure, all distributions converge to a point where R g ~ 9.5 
(fm) at vtlt — * nri p (px — > 0), indicating that the 'true 
source size' can be measured at small px, while at higher 
px the size parameter is found to be sensitive to space- 
momentum correlations. As it was mentioned above, the 
deuterons are sensitive to the space-time-momentum cor- 
relations in the same fashion as the two particle correla- 
tion function. Indeed, these effects were also found in the 
study of two pion correlation functions p9|j30[ |. 

Randomizing (or aligning) the nucleon space and mo- 
mentum vectors changes the space-momentum correla- 
tions. It affects both the deuteron yields and momentum 
distributions pi] . This effect is vividly displayed in fig- 
ure 2, where proton (circles) and deuteron (squares) av- 
erage transverse momenta {px) are shown as a function 
of rapidity. Plots (a), (b), and (c) are for normal, ran- 
domized, and aligned cases, respectively; plot (d) is the 
results from a calculation without rescatterings among 
baryons (rescattering here means interaction with pro- 
duced particles). The solid- and dash- lines in (a) and (d) 
represent the values of (pr) f° r kaons and pions, respec- 
tively. 

At mid-rapidity the values of the average transverse 
momentum (px) of pions, kaons, and protons are (px) = 
0.4, 0.6, and 0.85 GeV/c, respectively. The value for 
deuterons is about 1.4 GeV/c, see Fig. 2(a). These 
values are very similar to the results in References of 
The difference in mean px between proton and 
deuteron decreases from about 550 to 150 MeV/c as one 
moves away from mid-rapidity to y > 4. After random- 
ization, Fig. 2(b), the splitting between deuterons and 
protons in the mean transverse momentum becomes con- 
stant ~ 150 MeV/c. For the aligned case, the differ- 
ence changes from 580 to 280 MeV/c from mid-rapidity 
to y > 4. The change in transverse momentum at mid- 
rapidity is about 30 MeV/c, but at y > 4, the difference is 
about 100 MeV/c compared to the normal case. Similar 
to the random case, the results of calculations without 
baryon rescattering (Fig. 2(d)) show a constant differ- 
ence between deuteron and proton transverse momentum 
of about 150 MeV. 

Given the distributions shown in Fig. 1 and Fig. 2, 
one may discuss the physics in terms of collectivity. To 
proceed, we evaluate the average transverse velocities of 
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pions, kaons, and nucleons. In the following analysis, the 
collective velocity is defined as: 



</3t> = ( 



tttut 



(0.2) 



Figure 3 (a)-(d) show the velocities as a function of 
rapidity for different cases and Fig, 3(e) and (f) depict 
the mean number of nucleon collisions as a function of 
rapidity. 

Fig. 3 (a) shows that on average particles of different 
types are moving together with a similar velocity. These 
results indicate a certain amount of collectivity in Au + 
Au central collisions at RHIC energies.^ On the other 
hand, the RQMD calculations predict the charged pion 
(n + + 7r~) to nucleon (p + n) ratio is about 10 near mid- 
rapidity for the central Au + Au collisions at the RHIC 
energy. Levai and Miiller argue that in such baryon-poor 
region the equal magnitude of pion and nucleon flow ve- 
locities can be established only at an earlier deconfincd 
phase J3^]. Such deconfined phase, of course, was not in- 
cluded in the present calculation. The other interesting 
observation in Fig. 3 is that at a given rapidity, the val- 
ues of the velocity are closely correlated with the number 
of rescatterings: the larger the number of collisions, the 
higher the velocity (see Fig. 3(a) and (e)). In the case 
of no-rescattering calculation (see Fig. 3(f)), the num- 
ber of collisions is about 2.5 and the velocity is about 
zero over the whole rapidity range (Fig. 3(d)). Recall, 
that, in the no-rescattering case, the splitting in the av- 
erage transverse momentum according to mass vanishes, 
see Fig. 2(d). Once the space-momentum correlation is 
altered, the collectivity is destroyed (Fig. 3(c)). 

Note that RQMD predicts the averaged transverse col- 
lective velocity of mid-rapidity about 0.6 at the RHIC 
energy (Fig. 3(a)) while at the SPS (y/s m 20 AGeV) 
the value is about 0.4-0. 45c. 

In summary, using a microscopic transport model 
RQMD(v2.4) and a coalescence after-burner, we have 
studied the transverse momentum distributions of pi- 
ons, kaons, nucleons, and deuterons in different rapid- 
ity regions for central Au + Au ion collisions at \fs — 
200 AGeV . Employing the deuteron as a probe, we 
have demonstrated that a large number of rescatter- 
ings leads to the space-momentum correlation at freeze- 
out and is responsible for the decrease of the ratio of 
N 2 (proton) /N (deuteron) as a function of nvr- Should 
new physics occur at RHIC energy, a modification of 
the space-momentum structure will manifest itself in the 



2 The same model calculations, for Au + Au collisions at 
lower energies (few GeV per nucleon), indicate no such col- 
lective behavior although protons and light nuclear cluster 
distributions do show characteristics of collective motion at 
SIS and Bevalac energies (see |3l| and references therein). 



deuteron yields and its transverse momentum distribu- 
tions. These distributions can be measured in the STAR 
TPC and other RHIC experiments. 
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FIG. 2. Deuteron (filled square), nucleon (filled circle), 
kaon (solid-line), and pion (dashed-line) mean transverse mo- 
mentum as a function of rapidity. 
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FIG. 1. Gaussian radius parameter R g as a function of 
proton transverse mass bit for Au + Au central collisions 
(b < 3fm) at yfs — 200 AGeV. Filled circles represent the re- 
sults for original RQMD events. Open symbols represent the 
results for the events where the correlation between fir and 
r~r has been altered. Only mid-rapidity {\y\ < 1.0) nucleons 
are used for the plot. 
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FIG. 3. RQMD mean velocities (/3t) of pions, kaons and 
nucleons and mean number of collisions (ncl) for nucleons as 
a function of rapidity. 
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